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Abstract 



We consider the limiting location and limiting distribution of the largest eigenvalue in real 
symmetric (/3 = 1), Hermitian (/3 = 2), and Hermitian self-dual (/3 — 4) random matrix models 
with rank 1 external source. They are analyzed in a uniform way by a contour integral rep- 
resentation of the joint probability density function of eigenvalues. Assuming the "one-band" 
condition and certain regularities of the potential function, we obtain the limiting location of the 
largest eigenvalue when the nonzero eigenvalue of the external source matrix is not the critical 
value, and further obtain the limiting distribution of the largest eigenvalue when the nonzero 
eigenvalue of the external source matrix is greater than the critical value. When the nonzero 
eigenvalue of the external source matrix is less than or equal to the critical value, the limiting 
distribution of the largest eigenvalue will be analyzed in a subsequent paper. In this paper we 
also give a definition of the external source model for all /3 > 0. 

1 Introduction and statement of results 
1.1 Introduction 

In this paper we will be concerned with the distribution of the largest eigenvalue ^max{n-) in the 
following ensembles of matrices {M}: 

• The set of n x n real symmetric matrices, with the probability distribution function (p.d.f.) 



Pn,i{M)dM : 



g-nTr(yi(M) 





• The set of n x n Hermitian matrices, with the p.d.f. 



Pn,2{M)dM : 



1 



^-nTriViiM) 



(2) 
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• The set of 2n x 2n self-dual Hermitian matrices, with the p.d.f. 

Here for each /3 = 1,2,4, Cn^p is the normalization constant, Vp{x) (or Vp{x)) is a real-valued 
function which grows fast enough, A„^^ (or A„^^) is a fixed n x n real symmetric matrix, n x n 
Hermitian matrix and 2n x 2n self-dual Hermitian matrix respectively. The function V^(x) (or 
Va(x)) is called the potential function and A„^/3 (or A„^^) is called the external source. The rank 
of the external source A„^^ (or A„^^) is defined to be the number of nonzero eigenvalues of A„^^ 
if ^ = 1,2, or half of the nonzero eigenvalues of A„^^ if /3 = 4. These ensembles are called 
real symmetric, Hermitian and Hermitian self-dual random matrix models with external source 
respectively. Throughout this paper, we address the three types of ensembles as /3-external source 
ensembles with /3 = 1,2,4 respectively. 

Note that in these three external source models, the distributions of eigenvalues of M are 
unchanged if A„j3 (or A„^^) is changed into QAn^pQ~^ (or QA^^pQ'^), where Q is in the orthogonal 
group 0{n), unitary group U{n) and compact symplectic group Sp{n) for /3 = 1,2,4 respectively. 
Since we are only concerned with the distribution of eigenvalues of M, we assume A^^g (or A„^^) 
to be diagonal without loss of generality. To make our presentation uniform for all values of /3, we 
let V(x) be a fixed function and 

A„ := diag(ai, . . . ,a„) (4) 

be an n X n diagonal matrix. We assume that Va(x) (or Va(x)) and (or An,/3) are defined 
from V{x) and A„ such that 

Viix) = W{x), V2{x) = V4{x) = V{x), (5) 

A„,i = -A„, A„,2 = A„, (6) 

An,4 = diag(ai,ai,a2,a2, • • • ,an,a„) (7) 

Writing a 2n x 2n self-dual Hermitian matrix into 2x2 blocks ( ) , we can express it 



i,t=i 



as a quaternionic Hermitian matrix {qst)^t=i whose s,i-entry comes from the s,i-block by 

a + bi c + di 



,. , . , — a + hi + cj -\- dk. (8) 
c + ai a — bij 

In the quaternion form, the p.d.f. of the Hermitian self-dual external source model is 

PnAM)dM := ^e-"^T'^(^''(*^)-^"'4*^)dM, (9) 

Cn,4 

where 

V^ix) = 2V4(x) = 2V{x), A„,4 = 2A„ (10) 

with An defined in The A„^4 defined in pO|) corresponds to the A„^4 defined in d?]). In 

Appendix Rl we use the quaternion form p.d.f. Q of the 4-external source model to streamline the 
derivations for all /3. 
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In this paper we concentrate on the rank 1 case, i.e., 



A„ = diag(a,0,...,0). (11) 

n-l 

Random matrices are powerful tools to simulate Hamiltonians of complex systems. Different 
types of random matrices, namely the real symmetric (aka orthogonal), the Hermitian (aka unitary) 
and the Hermitian self-dual (aka symplectic) ensembles are used for physical systems with different 
properties of time-reversal invariance pi]. For random matrix models without external source in 
all the three types of ensembles, i.e., rank /3-external source models, the distribution of the largest 
eigenvalue has been studied extensively for all the three /?. If /3 = 2, for all real analytic potentials 
V2{x) under mild regularity conditions, the largest eigenvalue ^max(^) approaches e, the right-end 
point of the equilibrium measure of V2(x) (see ()16p below), with probability 1 as n — >• oo, and the 
limiting distribution is the GUE Tracy- Widom distribution. See e.g. [Hj and [12]. If /? = 1, for 
real analytic potentials Vi{x) satisfying the "one-band" condition (cf. Condition [2] in Subsection 
11.21 below) and mild regularity conditions, the largest eigenvalue with probability 1 approaches e 
as n — >• oo, and the limiting distribution is the GOE Tracy- Widom distribution. See [30j. If /? = 4, 
similar result can be obtained and the limiting distribution is the GSE Tracy- Widom distribution. 
See [12]. 

The random matrix model with external source was proposed by Brezin and Hikami [9], |10) 
to simulate complex systems with both random part and deterministic part. Although in all three 
types of random matrix ensembles the random matrix model with external source can be defined, 
due to technical reasons, only the Hermitian (/3 = 2) type has been studied for general potential 
functions. See e.g. [7J and references therein. 

In |3j, the Hermitian random matrix model with rank 1 external source was studied for all real 
analytic potentials V2{x) under mild regularity conditions. For convex potentials, the universality 
of phase transition was proved. Let V2{x) be defined by V{x) as in ([5]). In the rank 1 2-external 
source model, with probability 1, as n — )• oo 




W(n)^ r , ^ ' (12) 



where a is the unique nonzero eigenvalue of the external source A„^2 = A^, and xo{a) is a continuous 
increasing function in a G (^y(e),oo) such that xo{a) — e as a — )• ^^'(e), see (fT6|) and ([29|) . If 
a < ^V'{e), the limiting distribution of '^max('^) is the GUE Tracy- Widom distribution, and if 
a > ^y(e) the limiting distribution is Gaussian. For the double scaling a = ^V'{e) + the 
limiting distribution is the generalized Tracy- Widom distribution. If the potential is not convex, 
then new phenomena may occur. The "critical value" may be less than ^V'{e), and there may be 
"secondary critical values". The largest eigenvalue ^ma.x{n-) may converge to two or more points if 
a takes such values. The results were also obtained by Bertola, Buckingham, Lee and Pierce in [5] 
and [6j independently. 

For real symmetric and Hermitian self-dual matrix models with external source, known results 
are limited to special potentials. Let Vi{x) and Vi{x) be defined by V{x) as in the rank 
1 1-external source model with Gaussian potential {V{x) = on the real line) and Laguerre 
potential {V{x) = x — clog(x) on half of real line) are studied in e.g. [2], [28], [T7] and [IT]. The 
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limiting location of the largest eigenvalue is given by formula (jl2p . the same as in the corresponding 
rank 1 2-external source model, where a is the nonzero eigenvalue of A„ and A„^i is defined by 
([6]). If a > |y'(e), then the limiting distribution of ■^max("-) is Gaussian, with variance twice of 
that in the corresponding Hermitian (/3 = 2) external source model. If a < ^V'{e), the limiting 
distribution of ^max(?T') is the GOE Tracy- Widom distribution. The rank 1 4-external source model 
with Laguerre potential is studied in |32j . where the limiting location of the largest eigenvalue 
is found to be given by formula (jl2p . the same as in the corresponding rank 1 2-external source 
model, where a is the nonzero eigenvalue of A„ and A„^4 is defined by d?]). If a > ^y(e), then 
the limiting distribution of (,nia,x{n-) is Gaussian, with variance half of that in the corresponding 
Hermitian (/3 = 2) external source model. If a < ^V'{e), the limiting distribution of ^max('^) is the 
GSE Tracy- Widom distribution. In [32j the limiting distribution of ^max("') when a = ^^'(e) is 
also obtained. 

In the recent preprint |8j , Bloemendal and Virag obtained the limiting distribution of the largest 
eigenvalue ^maxin) when the potential is Gaussian or Laguerre, for all (3 and for all a. When o 
is at or near ^V'{e), they described the limiting distribution function of ^max{n) via the unique 
solution to a PDE. The recent preprint [25] by Mo indicates a new approach to study the limiting 
distribution of (,uiax{n-) in the rank 1 1-external source model with Laguerre potential when a is 
at or near ^V'{e), see also [26]. The contour integral formula in \25\ Theorem 1] is equivalent to 
that of Proposition 12.21 in this paper with (3 = 1 and Laguerre potential (cf. Remark lA.lj) . In [25] . 
Mo further simplified the integrand in the contour integral formula, (see [251 Theorem 3],) and he 
applied it in the asymptotic analysis in [26] to obtain a result similar to that in [8]. In this paper, 
we take a different approach to apply Proposition 12.21 in asymptotic analysis. The reader may also 
compare our paper with the paper [3] by Benaych-Georges and Nadakuditi, where they considered 
a different kind of low rank perturbations of large random matrices. 

In this paper, we consider the rank 1 /3-external source models with general potential V^(x) (or 
Vg(x)) which are defined by V{x). The conditions satisfied by V{x) will be given in Subsection 
11.21 We find that the "critical value" is independent of /3, and for all /3 = 1,2,4 find the limiting 
location of the largest eigenvalue Cmax(f^) when a, the nonzero eigenvalue of A„, is not equal to the 
critical value. When a is greater than the critical value, we also find the limiting distribution of 

Cmax(^) • 

Besides the asymptotic results summarized above, in Appendix [A] we also have an algebraic 
result: the definition of the /3-external source model with general /3 > 0. Here we note that the 
analytic method presented in this paper can be used to study the rank 1 /3-external source model 
with general /3. 

1.2 Assumptions on V{x) 

Throughout this paper, we assume four conditions on V{x), the function in ^ and (fTOj) . The first 
is 

Condition 1. V{x) is a polynomial of degree 21 and with positive leading coefficient. 

To state the other three conditions, we need a few definitions. Recall that if V{x) is regarded 
as a potential function on M itself, there is an equilibrium measure n associated to V{x) defined 
by a certain variational problem. See ()14p and (jl5p . and the general references |29] and [13] . Since 
V{x) is a polynomial, fi is supported on a finite union of intervals. In this paper we assume that 
V{x) satisfies the "one-band" condition: 
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Condition 2. The equilibrium measure /x associated to V{x) is supported on a single interval 
J = [61,62]. 

For the function V{x) satisfying Conditions [T] and [2], the equilibrium measure has the form 



d/i := ^{x)xjdx = —yj (62 - x){x - bi)h{x)xjdx, (13) 
zvr 

where XJ is the indicator function and h{x) is a polynomial of degree 21 — 2. The next condition 
assumed on V{x) is 

Condition 3. The function h{x) in the formula ()13p of the equilibrium measure fi of has 
only non-real zeros. 

The equilibrium measure dfj, = ^{x)dx is characterized by the conditions 

2 j log\x - s\^{s)ds - V{x) = £ for X G J, (14) 

2 j log\x - s\'^{s)ds -V{x) < e forxGM\J. (15) 

The last condition assumed on V{x) is 
Condition 4. The inequality (|15p is strict. 



Remark 1.1. Conditions [THS] are assumed to apply Proposition 14.11 in our paper, and they are not 
used anywhere else in this paper. If Proposition 14.11 can be proved under weaker conditions, e.g. 
the conditions assumed in [221 Theorem 1] 0, these conditions can be weakened accordingly. 

Remark 1.2. Functions V satisfying all Conditions [iHl] also satisfy the assumptions of V in [3l 
Formulas (6)-(8)]. Thus all the results in [3j on y can be applied in this paper. 

Remark 1.3. If V{x) is a convex polynomial with positive leading coefficient, V{x) satisfies Condi- 
tions [IHS] by ^214 Proposition 3.1], and it is straightforward to verify that V{x) satisfies Condition 

m 

1.3 Preliminary notations 

To state the results in this paper, we need a few more notations. We follow the notational convention 
in [3] to denote the right edge of the support of the equilibrium measure 

e := 62, the right edge of J = [61,62], the support of the equilibrium measure fi. (16) 

The so called g-function is defined by 

g(z) := j log{z - s)^{s)ds, for 2: G C \ (-oo,e). (17) 
For a G (0, ^^'(e)), define c(a) as the unique point in (e, 00) such that 

g'(c(a)) = / = a. (18) 



^Mariya Shcherbina informed the author that Proposition 14.11 can be proved under the consitions assumed in [221 
Theorem 1] through private communication. 
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The properties of g(x) used in this paper is summarized below (see [U Formula (30)]). 

g'(x) > 0, g"(x) < forxG(e,oo), 

V{e) + l V'{e) 
g(e) = , g(e) = — — , limg(x) = 0. 

2 Z a;— >oo 

For a > ^y(e), define c(a) := e. We also define two auxiliary functions 

G{z) = G(z; a) :=g(z) - V{z) + az, 
YL{z) = H(z; a) := - giz) + az + i, 

for z E C \ (— cxD, e). Condition [J] of V and the condition (jl4p imply that for any a 

1 £ 
G(e; a) = H(e; a) = -^V{e) + ae + -, 

G(x;a) < H(x;a) for x G (e,oo). 
The convexity of g(x) on (e, oo) yields that for u > c(a), 

U'iu;a)=a- [^>0. 

J u — X 

Define the set 

Av '■= {a G (0,oo) I there exists x G (c(a), oo)such that G(x;a) > H(c(a);a)}. 

It is proved in [3^ Lemma 1.2] that Ay is an open, semi-infinite interval. From Ay we define 

ac = aic{V) := inf ^y. 

It is also proved in [3l Lemma 1.2] that ac > 0. 
Let 

Gmax(a) := max G(x;a). 

xS:[c{a) ,oo) 

Define the discrete set 

J'v := {a G [ac,oo) |Gmax(a) is attained at more than one point}. 

For a > Re and a J^y, there is a unique xo(a) G (c(a), oo) such that (cf. [3| Lemma 1.3]) 

Gmax(a) = G(xo(a);a). 

For a > Sic and a G J'v, there are r > 2 and c(a) < xi(a) < X2(a) < • • • < Xr(a) such that 

Gmax(a) = G(xi(a); o) = • • • = G{xr{a); a). 

We define the set of secondary critical values as J'v \ {ad (cf. [3, Definition 1.3] ). 
Remark 1.4. If the potential V is convex for x > e, Jv = 0. See [3l Remark 1.2]. 
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1.4 Statement of main results 



Let V{x) be a function that satisfies Conditions [THll For any n and /3 = 1, 2, 4, let tlie n-dimensional 
^-external source models be defined by p.d.f.s ([T]), ([2]) and ^ respectively, with potentials 
(or Vg(x)) given by ([5]) and external sources A„^/3 (or A„^^) given by (jl]), ([6]) and (I7|. We assume 
that An has only one nonzero eigenvalue a, as in (jlip . In each n-dimensional /3-external source 
model, let ^max('T-) be the largest eigenvalue of the random matrix. The theorems below are stated 
uniformly for all /3-external source models (/3 = 1,2,4). In the case /3 = 1, we assume that the 
dimension n is even. For /3 = 1 and n is odd, the theorems below also hold, and we discuss it briefly 
in Appendix O First we show the limiting location of the largest eigenvalue. 

Theorem 1.1. The following hold for each fixed a as n —t- oo; 

(a) If < a < Sic, or a = S-c = ^V'{e) and ac ^ ^Ty, Cmax('^) — ^ e with probability 1. 

(b) If a > slc and Slc Jv, then ^max(^^) — XQ{a) with probability 1, where xo(a) is defined in 



(c) If a > Sic o,nd a £ Jy , then there exist r > 2 and xi{a), . . . , Xr{a) defined in ()30p . Un- 
der the assumption that G"{xj{a)) ^ for all j = l,...,r, then (,raax{n) converges to 
Xj{a) with nonzero probability Pj,/3(0) for j = l,...,r. Here Pj,p{0) are defined in (I168p 
and Y.]=iPjA^) = 1- 

Remark 1.5. If a < 0, Theorem 1 1 . ]|(a) | still holds, and the method of proof is similar to that in the 
< a < slc case. Since when a < there is no interesting phase transition phenomenon for the 
distribution of the largest eigenvalue (while there is a similar one of the smallest eigenvalue) and 
the proof is long and parallel to the a > case, we skip further discussions about the a < case. 

If o > ac, we have the limiting distribution of the largest eigenvalue. If a is not at or near 
secondary critical values, we have the following result that strengthens Theorem II. l](b)[ 

Theorem 1.2. The following hold for a > slc and a J'y as n —)• oo. 

(a) If G"(xo(a); a) / 0, then for anyT €R 

lim P I e^ax(n) < xo(a) + ^ = I = HT), (31) 

-fG"(xo(a);a)n/ 



where $(r) := j'^^e 2? d^ denotes the cumulative distribution function of standard 
normal distribution. 

(b) If G(^'^)(xo(a); a) / and G'^^\xQ{a)] a) = for j = 1, . . . ,2k — 1 where k > 1, then for any 
, p., . . , / -fG(^^)(xo(a);a)n \-^^^"^ \ f^^e-^^d^ 

If a > ac is at or near a secondary critical value, we have the following result that shows the 
double scaling case and strengthens Theorem I l.]]rc)[ 
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Theorem 1.3. Suppose that oq > ac and gq E J'y. Assume that G{x;aQ) attains its maximum 
at r > 2 points xi(ao) < ^2(00) < ••• < Xr{ao) in {c{ao), 00), and G"(xj(ao); ao) / for all 
i = 1,2, . . . ,r, then for 

a 

a = ao-\ — , (33) 
n 

where a is in a compact subset o/M, we have 

lim P I emax(n) < Xiiao) + ^ = = I >JPi,/3(a) I +Pi,fs{a)<!>{T), (34) 

^-^G"{xi{ao))nJ 




where pi^p{a) (i = 1,. . . ,r) are defined in (1168^ and Pj,/3(a) = 1- Furthermore, pr^/3{a) — )• 1 

as a 00 and pi^^(a) — )• 1 as a — )• —00. 

Remark 1.6. The phenomenon of Theorem 11.31 occurs for some quartic potential V that satisfies 
Conditions [THl For example, V{x) = 0.02093x^ - 0.16736x3 + 0.37448x2 + 0.11418x. 

In the case that a > slc, a S J^y ^md G"(xj(a); a) = at at least one maximizer Xj{a) of G(x; a) 
in (c(a), 00), we show hereafter an example when the number of maximizers of G(x; a) in (c(a), 00) 
is r = 2. The result for general case is similar. 

Theorem 1.4. Suppose that ao > Be and oq S Jv- Assume that G(x;ao) attains its maximum 
at two points xi(ao) < X2(ao) iciao),oo), with G"(xi(ao); ao) 7^ 0, G*^2'^)(x2(ao); ao) / and 
G(^)(x2(ao); ao) = for j = 1, . . . ,2k - 1. Then for 

a = ao - 9/3 h -, where := - — --^ — , x , (35) 

n n ^ P X2(ao) - xi(ao) 

and a is in a compact subset of M, we have 



n—^oo 



lim^P I W(n) < XI (ao) + , = =pi,/3(a)$(r), (36) 

^_|G"(xi(ao))n/ 

, p,, , ^ , / -fG(2^)(x2(ao);ao)n r^/^^^^ \ , , J-00 ^^'^x 

jirn^P I Uax(n) < X2(ao) + I ^^^yi I T I =Pi,/3(a) +P2,/3(a) jB^-^Z^3r^> 

(37) 

wherepi^p[a) andp2,/3{a) are defined in (jl83p . and pi /3(a)+p2,^(a) = 1- Furthermore, P2.i3{o!) — )• 1 
as a 00 anc? pi_/3(a) — )• 1 as a — )• —00. 

Remark 1.7. When /3 = 2, the probabilities pj^2{a) and pj^a) should agree with the p^i\X^) ™ [3 

Formula (52)] and the p'f\{a) in [3l Formula (63)] respectively. It is not obvious that they are the 
same, and we give the proof in Appendix iBl 

The limiting distribution of the largest eigenvalue when a < slc, as well as the limiting location 
of the largest eigenvalue when a is at or near ac < ^V'{e), will be analyzed in a subsequent paper. 

The paper is organized as follows. In Section [21 we calculate the limiting p.d.f. of the largest 
eigenvalue in the rank 1 /3-external source model as n — )• 00, based on Proposition 12. 1[ In Section 
[3l we prove Theorems 11.11 11.2111.31 and II. 4[ The proof of Proposition 12.11 is in Section [H 
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The starting point of the asymptotic analysis in this paper is Proposition 12. 2| the contour 
integral formula of the largest eigenvalue Cmax(f^)- Since its proof is combinatorial, we postpone it 
to Appendix [Al In this appendix we also propose a definition of the /3-external source model for 
any /3 > 0. In Appendix |B] we show that the results in this paper agree with those in ^ when 
/3 = 2. 



2 The p.d.f. of the largest eigenvalue 

In this section we compute f^^^^{n)^ the p.d.f. of the largest eigenvalue ^max(f^) in the n-dimensional 
/3-external source model with rank 1, as n — )■ oo. For /3 = 1, n is assumed to be even. We also 
assume that the only nonzero eigenvalue of the external source matrix A„ is a > 0. Recall that J 
is the support of the equilibrium measure fi associated to V{x), e is the right end of J, and c(a) 
is defined in Subsection 11.31 In this section we compute/estimate f^^^^{n)i^) fo^ real u. To be 
concrete, let e be a small enough positive constant. In Subsection l2.lt assuming that c(a) > e, we 
compute /^niax(n)(^) f'^r u £ [e + e, c(a) — e] up to a constant factor C„^^. In Subsection 12.21 we 
compute /^njax(n)(^) for ^ ^ ['^(^) +e, e + e~^] up to the constant factor Cn,i3- For u S (e + e~"^,oo), 
u G (— oo,e + e) and u G (c(a) — e, c(a) + e) in case c(a) > e, we give an estimate of /^n^ax{n)(^) 
Subsection 12.31 Note that throughout this section, u is always a real number. 

To facilitate the computation of /5max(n)(^)' define some notations. For any m, define the 
probability measure on M"^ 

dnm,t3ixi,...,xm) := ^ 1 A(xi , . . . , x^) j'^ TT e" ^dxi . . . (38) 

where Z^.^g is the normalization constant. Suppose F{x\, . . . ,Xm) is an integrable function with 
respect to the measure /U^,/? defined in ([38]) . define the expectation of F with respect to by 



'E'm,l3{F{xi, . . . ,Xm)) ■■= / F {xi , . . . , Xm)dfJ,m,l3{xi . . . Xm) ■ (39) 

For u G M and w G C \ {—oo,u), define the functions in u and w 

(x,) , (40) 
( ^ 

Zm,l3{u,w) ■.= 'Em,(3 I ^'m,/? (^^l , • • • , X^; M, w) | X(-oo,«) (^^j ) , (41) 

where 

^ e-^^'~'-'^)(u-x-f 
P„,/3(xi,...,x„;m,'u;) := TT , (42) 

{W - Xj)P/^ 

and we take the principal branch of {w—XjY/"^ for w G C\(— oo, Xj). For x < u and w G C\(— oo, u), 
define the function in x with parameter u and w 

p{x] u,w) := —V{x) + 2 log(ii — x) — \og\w — x\. (43) 
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We have 



If w = u we denote 



p{x] u) := p{x; u, u) = —V{x) + \og{u — x). 
Then we can state the technical tool in the asymptotic analysis of this section: 

Proposition 2.1. Let u > e and w & C \ {—oo,u]. 
(a) Suppose 



(44) 
(45) 



w = u -\ — , where z = s + it is in a compact subset of C \ (— oo, 0], 
n 



we have 



Zn-l,l3{u,W) = e 2 J u-=o Rp{u) 



exp 



73 



n I p{x; u)d^{x) 



(46) 



(47) 



where Rp{u) is defined in (|192p . 
(b) Suppose 



w = wq + i—i=, where wq > u and t is in a compact subset of M, 



we have 



Zn-l,l3{u,w) = e " ('"O-":) 



'0 "i?^(n,wo)exp 



/3 



n / p{x;u,WQ)dfi{x) 



(48) 



(1 + ^(1)), 
(49) 



where i?^(u, tuq) is defined in ()19ip . 
(c) Let e be a small positive constant. For all w such that dist(t(;, (— oo,e]) > e and \w\ < , 



\Zn-l,l3{u,w)\< Zn-l,l3{u,w) = exp 

where the factor 0(1) is bounded uniformly in w. 



n / p{x; u, w)dfi{x) 



0(1) 



(50) 



This proposition is a corollary of a theorem of Johansson [21' Theorem 2.4], and we put off its 
proof to Section HI 

For the asymptotic analysis in this section, we define four types of contours: Sf^ Ilf , Ff and 
/g , where x is a real parameter and si, S2 and s are positive parameters. We assume S2 > si/\/2 
for ^si,s2 ^^^^ allow s = oo in Ff. The contours IIJ and ^si,s2 ^^^^ used in Subsections 12.11 and 
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respectively. The contours and represent the local parts of IIJ and Sf^ around the 
point X respectively, which will turn out to be the saddle point in the asymptotic analysis. 



{w{t) I t G 



where 



w{t) 



X + e i t 
X + e 4 Si + i{t 



if < t < si, 

if Si < t < S2 + (1 



V2 



)Sl, 



{V2 - 1)S1 + (1 + i)s2 -t if t > S2 + (1 - 



[w{-t) 

{w{t) \t eR} where w{t) 



if t < 0. 



x + it if < t < s, 

x+{l + i)s-t if t > s, 
w{-t) if t < 0. 



(51) 



(52) 



={w{t) \t £ [-s,s]} where tt;(t) = + (it - |t|)/\/2. (53) 

Ig ={w{t) I t e [-S, s]} where tt;(t) = j; + it. (54) 

See Figures [H [21 [3] and m for these contours. For any real number r, we define 

K,s,ir) = {z G Sf^,., I > r}, U^.ir) = {z e | > r}. (55) 



The asymptotic analysis in this section is based on the contour integral representation of the 
p.d.f. of the largest eigenvalue ^max(?^): 

Proposition 2.2. Let i^max(?^) be the largest eigenvalue in the n- dimensional rank 1 P-external 
source model for (3 = 1,2,4, where the potential Vg(x) (or Vg(x)j is defined by ^ from V{x), and 
the external source matrix A„^^ (or P^n,p) is defined by ([6]) (or ^) from A„ in (llljl with a > 0. 
Then for any integer n if P = 2,A and for even integer n if f3 = 1, 

where Cn^p is a constant, Zn_i^p{u,w) is defined in (|30|) . and C is either Sf^ defined in ([5T]) or 
defined in (|52|) wit/i x > u. 

The proof of Proposition 12.21 is in Appendix Rl 
2.1 Computation of f(^^,,(n)iu) when m G [e + e, c(a) — e] 

Assuming e + e < u < c{a) — e, we use the contour integral formula (j56p of /g,„ax(™)(''^) take the 
contour C in (I56p as H^^^ , which is defined in ([52]) . Here vri is a large enough parameter such that 
the inequality ([63]) holds. For w S FI^^"-*, we parametrize it by w = w{t) as in ([52]) with x = c(a) 
and s = TTi. 
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a; + e 4 si' 



X + e 4. Si' 



X + is 



Figure 1: The contour SJ^ 



X + e 4. s 



X + is 



X 



S2 



57ri 

X + e i s 




Figure 2: The contour Ilf 

X + is 



X 



X — IS 



Figure 3: The contour Tf 



Figure 4: The contour 1^ 



Let L be a positive number. For sufficiently large n, the contour I'j^J^^ is part of H^i^ ■ By 

Proposition I2.l|(b)| for w = c{a) + itj^/n G ^'jfi^^i integrand in the contour integral of 
satisfies 



_ ^ I ,j!!l-^lv2 ^^n( fp(x-,u,c(a)W(x)+ac(a)) Rl3{u,c{a)) ^ ^^^^^^ ^^^^ 



Zn-l,l3iu,w) 

Using the asymptotic formula ([57|) . we have the result 
— 1^ Z^_,,(,{u,w)-^^-—dw = ^eHlpM-))^^^^^^ (58) 

ZTTl Jj^-^) (w — U)PI^ yjn 

L/ ^/n 

where 



Rp{u,c{a)) 



(59) 



and ei(L, n) is small if L and n are large. To be precise, for all e > 0, there is an Li > such that 
for all L > Li, |ei(L,n)| < e for n large enough. 

For w ma. bounded subset of 11^^"^ Proposition 12. life)] and estimate the integrand 

of the contour integral of 



Zn-l,p{u,w) 



62' 



62 



nP{u,w) 



\w — u\/^/' 



Ml), 



(60) 
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where 



P{u, w) 



x) + a^w, 



(61) 



and the 0(1) factor is uniformly bounded. Below we show that P{u,w{t)) decreases fast enough 
as t increases and t > 0. By the symmetry of the contour S^^^^ about the real axis, we see that 
P{u,w{t)) decreases fast enough as t decreases and t <0. 

If w{t) G n^^j"'* \ I'j^J^^ with t > TTi, i.e., w is in the ray from c(a) + ivri to — oo, we have 



dPiujwlt)) f c(a) — X + TTi — t 

— afi{x) — a. 



(62) 



dt J \w{t) — 

If TTi is large enough, we have that for all t £ [tti, oo) and all x €z J, there exists ci > such that 



c(a) — X + TTl — t ^ 1 



\w{t) — x\ 



c a — x 



Cl. 



Hence we have 



dP{u,w{t)) 



< 



d^{x) 



dt J c{a) — X 

For < t < TTi, hke ([62]) and ([MD, we have 

t 



Cl + a = — Cl. 



dP{u,w{t)) 
Jt 



\'w{t) — x\ 



;dlJL{x) < —C2t < 0, 



(63) 

(64) 

(65) 



where C2 is a positive constant depending on tti. 

Let Li be a large enough positive number such that —Li < c(a) and the inequality ()75p holds. 
By Proposition I2.1](c)| and inequalities ([Mil and (f65]l , 



1 



2TTi ./(n^(-)(-L,)U2';^)nc+ 



Zn-l,l3iu,w) 



€■2 



{w - uf/"^ 



dw 



1 



n 



exp 



'/3 



nP{u, w{L/y/n)) 



0(1), 



(66) 



where the factor 0(1) is bounded uniformly in L. Substituting the Taylor expansion ()210p of 
p{x;u,w) into ([66]) . we find 



,^nP{u,w{L / ^/n)) _ „§"■(/ p(x;u,c{a))d^i{x)+ac(a)) 



4 J 



(c(a)-xy2 _|_ 



(67) 



We write like IM 



27ri ./(n^(-)(_Zo\/^(;^)nc+ 



Zn-l,l3iu,'w)- 



dw = 

n 



§n(J p{x;u,c{a))dfj.{x)+ac{a)) ^^^j^^ 



By (j66p and (j67p . we find that e2{L,n) is small if L and n are large. To be precise, for any e > 0, 
there exists L2 > such that for all L > L2, |e2(L,n)| < e for sufficiently large n. 
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If w ^ (n^i'*^ \ n^j'^^(— Li)) n C+, i.e., w is in the ray from c(a) + ivri to —oo and < —Li, 



since — Aj| > tti for all Aj G M, by (j42p we have 



|P„_i,/3(xi, . . . , w)\<n, J] e" - xjf. (69) 

Hence substituting (|69l) into (|40l) . we have 



< ^ E„_i^;3(F„_i^^(xi,...,x„_i;'u)), (70) 

where 

F„_i,^(m, . . . ,x„_i;n) = r maxi<,<„x, < u, ^^^^ 

10 otherwise. 

Similar to (pDUj) and (pUTj) . we have (see Remark [H]) 

E„_1,;3(F„_1,^(X1, . . . , Xn-i;u)) = e§«/-V'(-)+21og(n-.)dMWQ(-^)^ (^2) 

^n-i,/3(w,u;) =cT2"^"ei"^-^(^)+2i°g("-^)'^^WO(l), (73) 

where in ([73|) ii; G (H^i"'' \n^i"^(— Li))nC+, and the 0(1) factor is bounded uniformly in w. Taking 
s = —^w, by substituting the estimate (f73]) of Zn-i^piu^w) into (i56]l . we have 

=e2'^(/~^(^)+21og(M-a;)rfAt(a;)-log(7ri)-aLi)Q^^-l-j 

where the 0(n^^) factor is uniform in Li. For large enough Li, 



-F(x) + 2 log(ti — x)(i/i(x) — aLi — log TTi < j p{x;u,c{a))dfi{x) + ac{a). (75) 
Substituting ([75|) into ([9^ . we have 



1 /■ p'^anw Is/ \ 

— / Z„_i,;3(n,«;)- -gj^dw = ^ei"(/ p(-;«.-W)'^M-)+-(a))^(i)^ (7g) 

The results of (|68l) and (|76]l give an estimate of the contour integral (|56l) on (II^i'^^ \/2^^'*^)nC+. 

The integral on (II^i'^^ \ I^J^) H C_ is similar, since both the contour and the integrand are 
symmetric about the real axis. 

Therefore, by ([58|) . (f68|) and (j76|) . we have that for u £ [e + e, c(a) — e], 

=C„^^_Le-§"/ V(x)<iMx)g§n(-V'{n)+2g(«)-^)gfnH(c(a);a) ^ (77) 
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where g(M) is defined in p!7|) . H(c(a);a) is defined in ([2T]) . and 



(78) 



2.2 Computation of f^^^^(n){u) when u E [c(a) + e, e + e ^] 

Assuming c(o) + e < u < e + e^^, we use the contour integral formula ([56]) of f^^^^{n)i^) take 
the contour C in ()56l) as SJ^^"^ ^ , which is defined in (I5ip . Here cii is a small enough parameter and 

(T2 is a large enough parameter, such that the inequalities ([55]) and ([ 

parametrize it by = as in (f5T]) with x = ii + n~^, si =(Ti and S2 = cr2 

Let L be a positive number. For sufficiently large n, the contour T^^" is part of S"^^^^ ^ • By 
Proposition 12 . ]|Fa)| foru; = u + 2;/nG T^l"^ the integrand in the contour integral of (|56|) satisfies 



hold. For w G 5^"+"^ \ "^^ 



(n,u;) 



62^ 



y J u — a; J 

(z/n)/3/2 



;i + o(l)). 



(79) 



Hence 



1 

27^ 



62' 



w — u 



)/9/2 



72' 



2/3/2 



-dz(l + o(l)). (80) 



Using Hankel's contour integral expression of Gamma function (See 6.1.4]), we find 



1 



djj,{x) 
u — x 



,B/2 



-dz 



- a 



dii{x) 



u — x 



(81) 



Comparing the integral on the right-hand side of ()80p with the left-hand side of ()8ip . we write 
analogous to ([58]) that 



1 



27ri 
where 



/ b(u, w) 



62' 



(tu — n)^/' 



-du; = n^/2-iei"(/f (^;")'^'^(^)+H7W^(n)(l + e3(L, n)), (82) 



- a 



u — x 



Til 



(83) 



and the term 63 (L, n) is small if L and n are large. To be precise, for any e > 0, there is an L3 > 
such that for all L > L3, |e3(L,n)| < e for sufficiently large n. 

Let L2 be a large enough positive number such that — L2 < u — a\l\pl and the inequality (j95|) 
holds. For w G S^+22"'(-^2) \ T^/^"', we use Proposition and find that ^ stih holds. 

Yi w = w (t) G and < t < o"!, i.e., w is in the line segment between u + n ^ and 



u + n + e i ai, we have 

dP(n,t£;(t)) 



(u-x-Fn-^)/\/2-t 
|t«(t) — 



dfi{x) 



a 

72- 



(84) 
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By (j24p . we know that H'(M;a) is a positive number. If cji is small enough, for all t G [0, cxi] and 
all X G J 

(85) 



{u- x + n-'^)/V2-t ^ 1 f 1 



U — X 



V2 



\w{t) - ^ \/2 V' 

Hence for < t < cJi, substituting ([85]) in ([8l|) . we find with the help of ([22 

dP{u,w{t)) ^ 1 fdf,{x) fl 1\ a 

= - -^H'iu; + - a) = -^H'(n; a). 

If w = w{t) ^ ^ s-iicl t > (72 + (1 — \/2/2)cri, i.e., u) is in the ray from u + n'^^ — (Ti/\/2 + i(T2 

to — oo, like ([62]) we have 



(86) 



dP{u, w{t)) _ [ u- x + n^^ + (T2- (\/2 - l)ai - t 



dt 



\w{t) — x\ 



dfi{x) — a. 



(87) 



If (T2 is large enough, like (f63]l we have that for all t G [172 + (1 — \/2/2)o-i, 00) and all 2; G J, if n is 
large enough 

n - X + n"^ + (T2 - (\/2 - l)(Ti - t 1 1^^,, , 

p-7- 12 < + T^H u; a . 

— u — X 2 

Substituting dMD into (HZD, we find that like ([88]), for t > 0-2 + (1 - \/2/2)(Ti 
dt <y ^ + 2^^^' 



a) — a = — -H'(u; a). 



(89) 



If = w{t) G and ai < t < (72 + (1 - V2/2)ai, i.e., is in the line segment between 

u + + e~cri and u + n""*^ — + i(72 we have 



dP{u,w{t)) f d 



dt 



dt 



log|w(t) — x\dfj,{x) < 0. 



(90) 



Thus by deO]), dMD, (HID and dlOD, we find that for a fixed L2, similar to jM 



1 



< 



(-L2)\r^+„" ^)nc+ 

2 /■Z+?i+n-l+o-2 + ('/2-l)cri 



^n-l,/3(M,^«)7 rgj^ 



■27r 



'L/n 



g|an«)(t) 

Z„_i,^('U, W^(t))^7T^ Tgrt 

{w[t) — UjP/^ 



2tt 



'L/n 



n exp 



—nP{u, w{L/n)) 



1 



1 



\w{t) - u\fil' 



;0{l)dt 



(91) 



=n 2 exp 



■/3 



nP(u, w{L/n)) 



w{L/n) - u\PI^ 
1 



0(1) 



|l + e^L|/3/2 



0(1), 
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where the last factor 0(1) is bounded uniformly in L. Substituting the Taylor expansion (j207p of 
p{x;u,w) into (f6T]) . we find that like (f671) 



^^nP{u,w{L/n)) ^ ^^n[J p{x;u)djA.{x)+au) | g" ^2 ^I^"72 + o(l) j . (92) 

Like (f68l) . we write 



(93) 

By (j9ip and (j92p . we find that e4(L, n) is small if L and n are large. To be precise, for any e > 0, 
there is an L4 > such that for all L > L4, |e4(L,n)| < e for sufficiently large n. 
Like (f74l) . we have 

1 /• gfanw) 

7r~ / 1 1 ^n-i,/?!^*,^^); W^*^^ 

Q^-n-H -V{x)+2\og{u-x)diJ.{x)~\og{<72)-aL2) Q(^^-l-^^ (g^^ 

where the 0(n~^) factor is uniform in L2. For large enough L2, 

-y(x) + 21og(ti — x)dfj,{x) — log(cr2) — 0X2 < J p{x;u)dfj,{x) + au. (95) 
Substituting ([95]) into ([Ml) , we obtain 

IT- . . Z„_i,^(n,u;)-^^-^dw; = n^/2-ie^(/P(^^")'^^(^)+Ho(l). (96) 



The results of ([93]) and ([96]) give an estimate of the contour integral in ([56]) on (E^+JJ^ ^ \ 
■pu+n ^) pi c The contour integral on (S"+'^~^ \ Lrt" ^) n C_ is similar, since both the contour 

L/n ' ^ a \ (71, (T2 \ L/n ' ' 

and the integrand in ([56p are symmetric about the real axis. 

Therefore, by ([MI), ([12]), ([M]) and ([96]), we have that for u £ [c{a) + e,e + r^]. 



/W(n)(") =C'„,^e-^^(")n^/2-ie^"(/ (A^^(n) + o(l)) 

=C„,;3e^«("^")(A^^(n) + o(l)), 

where G{u; a) is defined in ([20]) . 7W^(u) is defined in ([83]) and C„^/3 is defined in ([78 



2.3 Estimation of f^n-,^^{n){u) when u G (c(a) — e, c(a) + e), u > e + e^^ or u < e + e 
2.3.1 M G [c(a) — e, c(o) + e] or u G [e — e, e + e] 

In this subsubsection we use the inequality that if ui < U2 and w £ C \ (—00, U2), then 

Zn-l,l3iui,w) < Zn-l,ji{ui,w) < Zn-l,p{u2, w) . (98) 
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The inequality ([98]) is a straightforward consequence of the definitions (P0|) and (|1T]) of Zn-i,/3{u, w) 
and Zn_i^i3{u,w). 

For c(a) — e < u < c(a) + e, we use the contour integral formula of f^^^^{n)i''^) and take the 
contour C in (I56p as 11^2"^^^'^ i which is defined in ()52p . Here 7r2 is a large enough parameter such 



that the inequality puip holds. For w G II^j"^^'^'^, we parametrize it by w = w{t) as in (j52p with 
j; = c(a) + 2e and s = 7r2. From ([98]) we have for all u G [c(a) — e, c(a) + e] and L3 > — (c(a) + 2e) 



2vri ./n:W+2^'(-L3)nc+ 



62' 



< 



{w - n)/3/2 



27r 



^n-i,/3(c(a) + e,u;(t)) 



62 



f an5R(«)(t)) 



-dt (99) 



For w G 11^2"^'''^'^, we use Proposition I2.]|rc)] and find like ([6OI 



^n-i,/3(c(«) + e,?w(i))- 



,|anK(«;(t)) ^ 



62 



nP(c(a)+f,io(t))_ 



0(1) 



(100) 



Like (|64p . we have that for 112 large enough, for all t > 1^2 and x G J, there exists e'^^ > such that 
(cf. (IMl) and (IMD) 

dP{c{a) + e,?i;(t)) 



< -c'l. 



For all < t < vTi and x G J, we have like (j65p and (|90p that 

dP{c{a) + e,i(;(t)) 



< 0. 



(101) 



(102) 



Let L3 be a large enough positive number such that the inequality (llOSp holds. By (llOOp . (llOip 
and (fT02|) . we find hke ([MI) 



2 /•L3+c(a)+2f+7r2 _^ g|an5R(«;(t)) 

Z„_i,/3(c(a) + e,w{t))- — — TTTTKiit 



2ir 



w{t)-u\P/'^' 

n(J p{x;c{a)+l,c{a)+2l)d^i{x)+a{c{a)+2i)) Q^-^y (^]^Q3) 



62' 



Like ([Tip and (fMp . we also have 



27ri ./(n^(^a)+2f^(ncw+2f(_^^))^^_^ 



Zn-i,i3iu,w) 



a 



62 



{W - n)/3/2 

(/-V'(a;)+21og(«-x)d/i(x)-log(7r2)-aL3)Q^^-l^^ (104) 



where the 0{n ^) factor is uniform in L3. For large enough L3, 
y - y(x) + 2 log (n-x)(i/i(x)- log 7r2 - aLs < j p{x;c{a) + e,c{a) + 2e)dfi{x) + a{c{a) + 2e). (105) 
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Substituting (jlOSp into (|104p . we have 



I Zn-i r{u,w)- —j-dw 



"2 \\'-'ir2 

e2' 



.(/p(x;c(a)+?,c(a)+2e)c(At(a;)+a{c(a)+2?))^^-|^-j^ (-j^Qg^j 



The results of (|103p and (jl06p give an estimate of the contour integral ([56]) on Il^lf^^^'^ C+. 
The integral on n^j"^^^*^ ^ C!_ is similar, since both the contour and the integrand are symmetric 
about the real axis. Thus we obtain 



= C„^nl-'^/2g§nH(c(a)+2?;a)gfn((-y(c{a)+?)+2g(c(a)+e-)-^)+{y(c(a)+e-)-y{«)))Q(^j) 



For e — e<n<e + ewe use the contour integral formula (f56]) of /jn,ax(n)(^) i^ske the 

contour C in (f56l) as Xtl^^'^^'^'^ ^ the same as in the c(a) — e < n < c(a) + e case. We also apply the 
inequality (j98|) , and like pOSp have for all u G [e — e , e + e] 



1 



< 



Y /■c(a)+2f+L3+7r2 g|an5R(«)(i)) 

;r / Z„_i,^(e + e,^x;(t))— -y^dt. (108) 



Then we can find estimates similar to (|lU3p and pU6p . and obtain the estimate of f^^^^(n)i'^) similar 
to ()107p . We only state the result that for u G [e — e, e + e] 

/«.ax(n)M = Cn,/3nl-^/2ef"HW'^)+2.>)g|n((-y(e+?)+2g{e+e)-£)+(V'{^^^^^ (^gg) 

and skip details. 

2.3.2 n < e - e or u > e + e"^ 

In this subsubsection we first consider f^,^^^^{n)iu) for u < e — e. We use the contour integral formula 
(|56l) of /gn,ax(n)(^) ^^"i ^^^^ contour C as H^^'^. Here 7r4 is a large enough parameter such that 
the inequality (jll4p holds. 

We let Ci be any positive number, and define 

Cy = max(-y(x) + log(e-x)). (110) 



Let fi{x) be a function on M such that 



(1) fi{x) satisfies conditions (i) -(iii) mentioned in Proposition 14.11 

(2) On (-oo,e-e] 

h{^) = Cv (HI) 

and fe{x) is decreasing on (e — e, oo). 
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(3) 

J h{x)d^i{x) < -C (112) 



For w in the line segment of 11^+"^ from e + e — to e + e + i7r4, for aW x <u 

p{x; u,w) = — V{x) + 2 log(n — x) — log\w — x\ 

< - V{x) + 21og(e - x) - log(e - x) < Cy- 



(113) 



On the other hand, we assume that 714 is large enough such that for all x < e, 

- 1/(x) + 21og(e-x) -log^4 < Cy. (114) 

By (jllip and (|114p . it is straightforward to check that for all w in the two rays of n^^i^'^, from 
e + e + i'Ki to —00 and from —00 to e + e — 11:^ respectively, and for all x < u 

p{x;u,w)<Cv- (115) 
Then by (fTT3]l . (fTT5]l . (HQ]) and (|M]), we have for all w G n^+'' 

By Proposition O and (fTT2]) . 

En-i,i3 (^6 2 ^i=i = e-^'^^^-O(l). (117) 

Using the estimate ()116p and (]117p of Zn-i,i3{u,w), we find by direct calculation 



1 " 



e,^n(a{^+'^)-C~,)0{l). (118) 



Thus by ((Ml) 

/?_(n)(n) =C'„,^e-^^(")e^W^+^')-^^-)0(l) 

=C„/3n^-^/2e2"("^^")+"(''+')"^«+-^^('')'^''^''))0(l). 



(119) 



Note that can be any number, and the last 0(1) factor in (|119p is bounded uniformly for all 
u < e — e. 

Next we consider f\^^^(n){u) for n > e + e^^. We use the contour integral formula (f56]l of 
/^max(n)(^) and take the contour C as 11^+^, which is defined in ([5^ . 
Let 

Kim = miny(x), (120) 

and denote 

Cv ■■= J (-Klin - 2x)dfiix). (121) 
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For all w G n""*"^ and x < u, 



p{x] u,w) = — V{x) + 2 \og{u — x) — \og\w — x\ 
<-V{x) + 2\og{u-x) 
<-V{x) + 2{u- x) 



2x + 2u. 



Thus similar to (fTTBI) . we have by (fT22]) that for w G 11"+^ 
By Proposition 14.11 



0(1). 



Using the estimates (I123P and (11241) of Zn-i,j3{u,w), we have like (lllSp 



Thus by ([56 



1 /" gfnaw 
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n((2+a)M+a+C'v) 



(122) 



(123) 



(124) 



(125) 



(126) 



Note that the last 0(1) factor in (|126|) is bounded uniformly for all ti > e + e 

3 Proofs of Theorems [TI, [Q D and D 

In this section we prove the main theorems in this paper. We divide the proofs into three subsec- 
tions. In Subsection 13. H we consider the case that < a < ac and the case that a = ac = ^^'(e) 
and ac ^ ^-^d prove Theorem II. I] ( a) [ In Subsection 13.21 we consider the case that a > ac and 
a JV) and prove Theorems 11.21 and ll.l|(b)[ In Subsection 13.31 we consider the case that a > ac 
and a S Jv and prove Theorems II. 3| ll.l|(c)| and 11.41 



3.1 Proof of Theorem ll.]Ji(a) when < a < Hc, or a = Hc = \V'{e) and He ^ Jv 



First we consider the case that < a < ac. Let e be a small positive number, such that e + e < c(a) 
and e + e~^ > c(o). Furthermore we assume that e is small enough such that the inequalities (jl27p . 
(fT29D and (fT3TD hold. 

The condition < a < ac implies the inequality H(c(a);a) > Gniax(a), see ([25]) - ([30]) and O 
Lemma 1.2(d)]. We assume that 



H(c(o); a) 



c(a) > e, 



(127) 
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where Gmax(a) is defined in ([27]) . Then by ([97|) with e = e, for all u G [c(a) + e, e + e ^] we have 

/C.ax(n) W = C„,;3ei"^— ('^)O(l) = C„,;3e^""('=('^)''^)0(e-i^"). (128) 
We assume that for u,v £ [c(a) — e, c(a) + e] 

II{v + e; a) + (-^('u) + 2giv) - £) < H(c(a); a) - e. (129) 
Then by ()107p with e = e, for n G [c(a) — e, c(a) + e] 

/e_(n)(^) = C„,/3e^««'^)^'^)0(e-§-). (130) 
We assume that for u > e + e^^, 

- V{u) + {2 + a)u + a + Cv + j V{x)d^i{x) < H(c(a); a) + (e + er^ - u) - 1. (131) 

Then by (|126p with e = e, for n > e + we have uniformly in u 

/5_(„)(.x) = C„,;3e^"W'^)^'^)0(ei'^(^+^-'-«-i)). (132) 
Let C be large enough such that 

a{e + e)-C- V{u) + ^ V{x)dfi{x) < H(c(a); a) + e - e - u - 1. (133) 
Then by (I119P with e = e and = C, for u < e — e we have uniformly in ti 

/W{n)(n) = C„,^ef'^H(c(a);a)o(gfn(e-.-«-l))^ (^34) 

By ([MD, (fT30D . (fT32]l . ([IMD, we find that 

P(W(n) > c(a) - e or emax(n) < e - e) = C,,^e^'^HW'^);'^)0(e-i^"). (135) 
Let e' < e/2 be a small positive number such that 

• —V{u) + 2g{u) — £ is decreasing on [e, e + 2e']. 

• —V{u) + 2g(u) — £ attains its maximum on [e + 2e', c(a)] at e + 2e'. 

• Let e" := -{-V{e + 2e') + 2g(e + 2e') - £). Then e" < e. 
Then by ()77p with e = e', we have that 

P(emax(n) G [e + 2e',c(a) - e]) = C„,;3et"HW'^)''^)0(ef ^""), (136) 

and as n — )• cx) 

IP(gmax(n) £ [e + e^e + 2eq) ^ 
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The probabilities (|135p . (jl36p and (jl37p imply that the conditional probability 

P( G [e + e', e + 2e'] | U^{n) G M \ (e - e, e + e')) ^ 1. (138) 

Since [e + e', e + e'] U (e - e, e + 2e') G [e - e, e + e], (fT38|) implies that 

H^mM G [e - e, e + e]) ^ 1. (139) 
Taking e arbitrarily small, we prove Theorem ll.]|(a)| when < a < a^. 

The case when a = ac = ^^'(s) and ac ^ is similar. Let e be a small enough positive 
number. Since a = ^V'{e), we have c(a) = e. Since a = ac ^ c7y, there exists e' > depending on 
e such that for all u > e + e, G(n; a) < H(e; a) — e'. Thus like (I128[) . for u G [c(o) + e, e + e~^] we 
have 

/e_(„)(u) = C„,;3e^"(''^")0(e-§^'"). (140) 

When e is small enough, (jl32p and (|134p also hold. Then by arguments similar to (|135p - (|139p . we 
prove Theorem II . ]](a) | when a = slc = ^V'{e) and ac ^ 



3.2 Proof of Theorems 11.21 and 11.1 1(b) when a > ac and a ^ Jy 

Let e be a small enough positive constant such that the maximizer 2;o(a) of G(x;a) in [c(a), oo) is 
less than e + e~^, and the inequalities (|14ip . (jl45p and (|147p are satisfied. 

First we consider the case that ac < a < iy(e), i.e., c(a) > e. We assume that 

Gmax(a) -H(c(a);a) > e. (141) 

Then by (i77|) with e = e, for G (e + e, c(a) — e) we have 

/€.axW(") = C„,;3n(i-«/2el"HWa);a)Q(^) ^ C„,^ef 0(6" ^") . (142) 

We assume that for all n G [e — e, e + e] 

H(c(a) + 2e; a) + (-^(e + e) + 2g(e + e) - ^) + (^(e + e) - V(u)) < Gniax(a) - e. (143) 

Then by (jlOOp with e = e, for u G [e — e, e + e] we have 

4„ax(n)(^) = C„,;3ef"«— ('^)0(e-§-). (144) 
We assume that for all u G [c(a) — e, c(a) + e] 

H(c(a) + 2e; a) + (-y(c(a) + e) + 2g(c(a) + e) - £) + (T/(c(a) + e) - ^(u)) < G^^ax - e. (145) 
Then by (|107p with e = e, for u G [c(a) — e, c(a) + e] we have 

fu.An)i^) = C„,^e^G— ('^)0(e-§-). (146) 
We assume that for n > e + e^^, 

- Viu) + {2 + a)u + a + Cv + j V{x)dfi{x) < Gmax(a) + (e + e'^ -u)-l. (147) 
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Then by ()126p with e = e, for n > e + e ^ we have uniformly in u 

4_(„)(u) = C„,^ei"«— (")0(ei"(^+^''-'^-i)). (148) 
Let C be large enough such that 

a{e + e)-C- V{u) + j V{x)dfi{x) < Gmax(a) + e-e-u-l. (149) 

Then by (I119P with e = e and = C, for u < e — e we have uniformly in u 

/^_(„)(^) = C„,^ei'^^—(")0(et"(^— «-!)). (150) 
By (HUD, (HMD, (fTiSl) and (fTSpp . we find that 

P(emax(n) > e + or emax(n) < c(a) + e) = C„,^e^^— (")0(e-i^"). (151) 



In the case that a > ^y(e), i.e., c(a) = e, we find that inequalities (I144p . (I148P and (|150p still 
hold, and the estimate ()15ip holds with c(a) = e. 

For M E [c(a) + e, e+e~^], we have the asymptotic formula (f97D and xo(a) is the unique maximum 
of G{x; a) in [c(a)+e, e+e^^]. If we further assume that G"{xQ{a); a) / 0, by the standard Laplace's 
method we have that 



e [c(a) + e,e + e-i]) = .J"",, ■M^(xo(a))C„,;3ei"^— ^ (1 + o(l)), (152) 

— ^G"(a;o(a); a)n 



2 



and for any T £ 



^max(?^) e 



c(a) + e,xo(a) + 



fG"(xo(a);a)n 



p^ J"", , , ■M;,(xo(«))C„.,^e^^— W(^(r) + o(l)). (153) 
-^G"(a;o(a); a)n 

The probabilities p^Tp . (fT52D and (fT53D imply Theorem [L^a)] 

If the second derivative of G{x;a) vanishes at xo(o), due to the analyticity of G{x;a), there 
exists k > 1 such that G^^\xo{a);a) = for j = l,...,2fe - 1 and G^^''\xo{a);a) / 0. By the 
Laplace's method we have 



P(W(n) G [cia) + €,e + €-']) 



{2k)\ 



l/{2k) 



fG(2^)(xo(a);a: 



2 

f>00 



X / e-"''dxA1;3(xo(a))C„,;3et"f^— ('^)(l + o(l)), (154) 



— oo 
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and for any T € M 

-l/(2fe) 



c(a)+6,xo+( ^^y^ 1 T 



\ lj{2k) „ 



-fG(2fc)(xo(a);a) 



n 



e-^''dxA1/3(xo(a))C„,;3e^"^--('*)(l + o(l)). (155) 



The probabilities (I15ip . (I154p and (1155P imply Theorem ll.^(b)[ Finally, Theorem 11.21 implies 
Theorem ll.l](b)[ 



3.3 Proofs of Theorem [HH [njl(c)] and [TH when a > and a G Jy 

Let ao > ac, and e be a small positive constant such that the inequalities ()142p . (|144|) . (I146|) . (I148p 
and ()150p hold with a = oq- It is easy to verify that there exists a positive number e depending on 
e such that if we take a = oq + e' with e' G [-e, e], the inequalities (I142p . (I144p . (|146p . (11481) and 
()150p still hold with the same e. Thus the estimate of probability ()15ip still holds with a = oq + e'- 
If we further assume that 

c(a) >c(ao) - |, (156) 

Gmax(a) >Gmax(ao) " |, (157) 

by ()15ip we obtain that 

nUUn) > e + or U-M < c{ao) + |) = C„,^e^^— ('^«)0(e-^"). (158) 

First we assume that oq G ■Ty and G(x; oq) has r > 2 maximizers xi(ao) < ^2(00) < • • • < Xr{ao) 
in (c(ao), 00), and all of them are less than e + e~^. Further we assume that for alH = 1, . . . , r 

G"{xi{ao)) + 0. (159) 

We take 

a = ao + -, (160) 

n 

where a is in a compact subset of M. Since a G [oq — e, oq + e] for sufficiently large n, the estimate 
(I158P is applicable to a. 

For X around Xi(ao) (i = 1, . . . , r), we denote 

X = Xi(ao) H , (161) 

y-f G"(xi(ao);ao)n 

For in a compact subset of M and x given by (|16ip . we have 
-G(x; a) =^G(x,(a); a) - ^- + 0(n-2) 



2 ' ' ' 2 ^ ' 2 n 

= -G'max(ao) + -axi(oo) — 

2 2 n 2 n 



^^^^^ 

77Gmax(ao) + ^axj(oo)- - ^- + 0(n"^). 
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Let ej (i = 1, . . . , r) be small enough constant numbers such that Xi(ao) is the unique maximum of 
G(x;ao) in [xi(ao) — ei,Xi{ao) + ej]. Applying the standard Laplace's method to (1971) . near Xj(ao) 
(i = 1,2), we obtain that 



'(Cmax(n) e [xiiao) - €i, Xi{ao) + ei]) = e 



^2 



^ J^;— ^f— — -M/3(^.(«o))C„,^e§'^«— ('^o)(l + o(l)), (163) 
y -f G"(xi(ao);ao)n 



and for any T in a compact subset of M, 

T 



P emax(n) e 



Xi{ao) - ei,Xi{ao) + 



^ -^G"{xi{aQ);aQ)n 
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§...(ao)+^ / y A^^(x,(ao))C„,;3e^^— ^"°)(c|>(r) + o(l)), (164) 

y -f G"(xi(ao);ao)n 

There exists e" > depending on ei, . . . , e,- such that for sufficiently large n 

G{x; a) < Gmax(ao) - e" for x G [c(ao) + |, e + e"^] but x (2;i(ao) - e^, a:i(ao) + e^). (165) 

Then we find that the probability 

P(^emax(n) e [c(ao) + |,e + e-i]\|J(x,(ao)-e,,x,(ao) + eO^ = C„,^e^«— ('^«)0(e-§^""). 

(166) 

For i = 1, . . . , r, let 



B,^p{a) := e§°"'("°)J ^ :A^/3(x.(ao)) (167) 

and 



-f G"(3;i(ao);ao) 



PiA^) ■= v^? „ / N - (168) 

From (jl68p we immediately find Yl^j=i Pj,l3{^) — 1- By (|167p and (jl68p we find that limQ_5.oo Pr,/3(o) = 

1 and limQ,^_ooPi,/3(a) = 1. 

The probabilities (fT58]) . pM|) and show that for any T € M and i = 1, . . . , r 

W(n) < Xiiao) + , ^ I = I J]Pi,/3(«) I +PiA(^){HT) + o{l)). (169) 

^-fG"(x,(ao);ao)n/ \i=i / 

Therefore Theorem 11.31 is proved. Theorem 1 1 . ]](c) | is a consequence of Theorem 11.31 with a = 0. 

Next we consider the case that r = 2 and for k > 1 

G"(xi(ao)) /O, (170) 
G(2'^)(x2(ao)) /O, G(-')(x2(ao)) = j = 1, . . . , 2A; - 1. (171) 
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and take (see (j35]) ) 



logn a , 2 i — ,^ 
a = Go- 9/3 h-, where qp ■= ^ — -, — r — r, (172) 



where a is in a compact subset of R. For x around x\{a)^ we denote 



x = xi(ao)H , , (173) 

|G"(xi(ao);ao)'^ 



and for x around X2(a), we denote 



/ -fG(^^)(x.(ao);ao)n r^^^^^^^ 
x = ^2(ao)+( ^^y^ I 6. (174) 



For ^1 in a compact subset of M and x given by (jl73p , we have hke ()162p 

^G(x; a) =^G(xi(ao); a) - ^- + o(n"^) 
2 2 In 

=T7Gmax(ao) + -a;i(ao) — -Tr- + o{n ) 

2 2 n Z n 

For ^2 in a compact subset of M and x given by (jl74p , we have 

^G(x; a) =^G(x2(ao); a) - if- + o(n-i) 
2 2 n 

2 2 n n 



(175) 



(176) 



Let ei {i = 1,2) be smah enough constant numbers such that Xj(ao) is the unique maximum of 
G(x; oo) in [xj(ao) — e^, Xj(ao) + Cj]. Applying the standard Laplace's method to (f97ll . near xi(ao), 
we obtain similar to (jl63p and ()164p that 

-|a:i(ao)g,3log?i 

P(W(ri) G [xi(ao)-ei,xi(ao) + ei]) = ^ e^^'^^"")" 



Y -|G"(xi(ao);ao) 
and for any T in a compact subset of R, 



A^^(xi(ao))C„,^e^^— ('^°)(l + o(l)), (177) 



^max(?^) G 



T 

xi(ao) - ei,xi(ao) + 



^J-^G"{xi{ao)■,ao)n 



^ et-^("°)" / ■M^(x2(ao))C„,^ef"^— (-°)(^(T) + o(l)). (178) 

V -2G"(a;2(ao);ao) 
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Applying the Laplace's method to (j97|) . near X2{ao), we obtam 

-■§X2(ao)g/3log7i 

G [x2(ao) - e2,X2(ao) + £2]) = e2-2H)a 



\ 1/(2A:) 
{2k)l \ 



-fG(2^0(x2(ao);ao) 



2 

and for any T in a compact subset of M 



A1^(x2(ao))C„,^e§"«— (^y" e^^'d^ + , (179) 



^max(ra) G 



X2(ao) - e2,X2(a) + I ^^^^y^ I T 



X 



/ e«2'(i6 + o(l)y (180) 
J —00 / 



Also there exists e" > depending on ei and e2 such that the estimate ()166p holds. The probabilities 
(HMD, (HMD, (fTTTl) . (fT75]) . i^M!), (fTSO]) show that the probability that Cmax(f^) is in [xi - ei, + ei] 
or [x2 — £2,^2 + €2] approaches 1 as n — )• 00. 
Let 



D,^f, :=e^Mao)» M^{x,{a,)), (181) 

y -f G"(xi(ao);ao) 

I)2,/3 :=e^^('^'')- . , A^/3(x2(ao)) H e^" dC2, (182) 



fG(2'=)(x2(ao);ao) 



and for i = 1,2 



L'i,/3(a) +L>2,/3(a) 

From (fTSSD we immediately find that Vx^ia) +P2,/3(a) = 1- By (fW]) . (fTH2]) and (HMD, we find that 
limQ,^ooP2,/3(a) = 1 and limQ,^_oo Pi,/3(a) = 1. 
Because 

We further find from (fTTT]) . (fTTHD . (fT79l) . (fT80]l that for any T G M 



emax(n) < xi(ao) + , ^ ^= I = Pi,^(a)($(r) + o(l)), (185) 



-fG"(xi(ao);ao)n 
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Cmax(ra) < X2(a) + 



-fG(2fc)(xo(a);a)n\"'^^''\ 



+P2,/3(a) ( j5l^^Z^2r^ + o(l) 1 • (186) 



Therefore Theorem 11.41 is proved. 

4 Proof of Proposition 12.11 

The proof of Proposition 12 . 1 1 is based on a theorem of Johansson [21] . For the convenience of readers 
we state it bellow. 

Proposition 4.1. !121\ Theorem 2.4] Suppose the function V satisfies Conditions and is 



the equilibrium measure associated to V. Let f be a real function that satisfies conditions (i)- (Hi) 
below, with s = 2 if (3 = 2 and s = 17/2 if P / 2. Then there are a quadratic functional A on f 
and a signed measure v on supp(^) = J = [61,62] which do not depend on n, such that as n —)• 00 

logE„_i,^(e^"=i ^(^^■)) = {n-l) j fix)df,{x) + (| " l) / fi^)M^) + ^Mf) + 0(1). (187) 
The quadratic functional A is defined by 

A{f) = ]^jj{x)5f{x)dx (188) 

where 5^ is given by 

^/(,) . '1 / / (189) 



2irV(i--6l)(62-i) ']., 

From the quadratic functional A, we define the inner product (•, ■)a by 

{f,9)A := \{A{f + g)- A(f) - A{g)). (190) 

The explicit formula of I'ix) is more complicated and is given in (21', Formula (3.54)]. The conditions 
mentioned in Proposition 14.11 are (see [21], Page 157]) 

(i) fix) < C{V{x) + 1) for some constant C, all x £R. 

(ii) |/'(x)| < q{x) for some polynomial q{x) and all x G M. 

(iii) For any xq > 0, there is an a > such that hip^o ^ H'^'^", where H^, s > 0, is the standard 

Sobolev space, and ipx^ G C°° is the function such that tpxoi^) = 1 if < 2^0, tpxoi^) = 
if > xq + 1 and < ipxoi^) ^ 1- 
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The function Rf3(u,w) appearing in Proposition 12.11 is defined by 



Rp{u, w) := exp 



^ ( ( ^ - 1 ) / p{x; u, w)dv{x) - I p{x; u, w)dfi{x) + A{p{x; u, w)) 



for u £ (e, oo) and w & C \ (—00, u). li w = u, we denote 

Rl3{u) := Rp{u,u). 
To facilitate the proof of Proposition 12.1^ we define 



(191) 



(192) 



,Xm;U,W) 



3! ' 



(193) 
(194) 



where c is a parameter no greater than u. When c = u, Zra,i3{u,w; c) and Zm^p{u,w;c) become 
Zm,p{^,w) and Zm,fi{u,w). 

The proof of Proposition 12.11 is as follows. Recall that e is the right edge of J, the support of 
the equilibrium measure. Let c = (e + u)/2. We write 



Zn-l,l3{u,w) = Zn-l,l3{u,W;c) 



Zn-l,p{u,W]c) Zn-l,p{u,w) - (u, W; c) 



Zn-i^0{u,w;c) 



+ 



Zn-i,i3{u,w;c) 



In case (a) where w is given by (I46p . we assume the results 

Zn-iAu, w) =e-^^ ^i?/3(u)ei"^P(^;")'^^(^)(l + o(l)), 
Zn-iA'^, w; c) =e-^ ^ '^i2^('u)ei"i'P(^;")'^^(^-)(l + o(l)), 

Zn-i,i3iu,w;c) 



Then we have 



Zn-l,l3{u, W) - Zn-l,l3{u, W, c) 



Zn-i,i3iu,w;c) 



< 



Zn-l,l3iu, w) - Zn^i^U, W, c) 
Zn-l,l3{u,W;c) 



(195) 



(196) 
(197) 
(198) 



0(1), (199) 



and find that (gT]) is the consequence of (fT96]l . (fWD . (fTMD and (fT99D . 
In case |(b)| where is given by (^8]) , we assume the results 

Zn-i,p{u,w) =e ~-' K^ii^(u,u;o)e2"/P(™»)'^'^(^)(l + o(l)), 

r diJ,{x) a 

Zn-i,p{u,w;c) =e~~^ (^^^^i?^('u, 'u;o)e^"i^f(™»)'^^(^) (1 + o(l)), 

Z„_i,^(ti,ty;c) -i&L^fM2± 

— — =e ^ •''^0-^+0(1). 

Zn-l,B{u.,W\c) 



(200) 
(201) 
(202) 
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We still have (fTMIl . and (09]) is the consequence of (l200D . (pOTTl . (12021) and (fT99|) . 

we assume 



In case (c) 



(203) 



Then we immediately obtain ()50p . 

Below we prove the asymptotic formulas (fTM]l . (fWll . (fTMD . (I200D . (l20T]l and (l203]l . 

Proo/ o/ (fT96D . (fWll . ([200]) . (f20TD and (I203D . In the proof, c stands for c = (e + n)/2 or n. 
For fixed u, w and c, let /i(x) and 72(2;) be two functions on M, such that 



(1) fi{x) and f2ix) satisfy the conditions (i) -(iii) mentioned in Proposition 14.11 

(2) fi{x) = f2{x) = p{x;u,w) for x < e. 

(3) fi{x) > p{x; u, w) for x G (e, c). 

(4) There exists xq G (e,c) such that f2ix) < p{x;u,w) for x G (e,xo) and fi{x) — f2{x) > log 2 
for X > Xq. 

As a consequence of the properties of fi{x) and /2(x), we have 



2E 



n-l,/: 



n-l 



<E„_i,/3 j \Pn~l,l3{xi, ■ ■ . ,Xn-i;U,w)\ X{-oo,c){^ 



(204) 



By Proposition 14.11 we have for both i = 1,2 that 



E„_i,^(e2^^=i^'(^^)) =exp 



/3 



(n — 1) / —p{x;u,w)dfi{x) + {— — l) / — p(a;; u, ty)(iz^(2;) 



/3 



2 . //3 , 



(1 + ^^(1)) 



/3 



n / p{x;u,w)dfi{x) 



-Rp{u, w) exp 



Thus by the sandwich inequality (j204p and (|194p we obtain 



/3 



(l + o(l)). 



n.-l 



(205) 



-^n-l,/3(tt,W^;c) =En-l,i3 ( |-P„_l,/3 (xi , . . . , X„_i ; U, U;) | JJx(-oo,c)(= 



(206) 



-Rf3{u, w) exp 



n / p(x; n, w)d^{x) 



(l + o(l)). 



By (I206p . we complete the proof of (j203p with c = u. Let w be given in (|46p . we have uniformly 
for all X < e that 

(207) 



p{x]u,w) = p{x]u) n ^ + 0{n ^). 



u — X 
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Thus 
and 



exp 



n / p{x; u, w)dfi{x) 



Rp{u,w) = Rp{u) + 0{n-^) 



_ /3s r d^{x) 

e "2" "-^ exp 



n / u)dfi{x) 



(208) 



(l + 0(n-i)), (209) 



and we obtain the proof of (jl97p with c = c and the proof of ()196p with c = u. Let w be given in 
(j48p . we have uniformly for all x < e that 



p{x; u, w) = p{x; u; wq 



Thus 

and 
exp 



2{wo — xY 
Rp{u, w) = Rp{u, wq) + 0(n"^) 



n-^ + 0{n-'^). 



(210) 



(211) 



/3 



n f p{x; u, w)dfi{x) 



■' (tllQ- 



exp 



/3 



n / p{x;u,wo)d^{x) 



(l + 0(n-i)), (212) 



and we obtain the proof of (j20ip with c = c and the proof of (|200p with c = u. 



□ 



Remark 4.1. By the same method, we can evaluate E„_i^^(F„_i^^(xi, . . . ,rc„_i;u)) where Fn-i,i3 
is defined in (1711). 



Proof of (jl98p and (j202p . We consider (—oo, c)" as a probability space with the probability mea- 
sure 

1 

|^'n-i,/3(a::i, • • • , a;n-i; ti, 'u;)|d^n-i,/3(xi, . . . , Xn-i), (213) 



Zn-i,i3{u,w;c) 

where dnn-i,i3{xi, . . . ,Xn-i) is defined in ([38]) . Let ^^.i^ be a random variable on (— oo,c 
such that 



,n-l 



(214) 



where the range of the argument is taken to be (— vr, vr]. 
We define 

1 ^ "^/^ 

o-„ := ^ ' 



V — X 



for any f > e. For w given in ()46p . we will show 

E(5',^-i,/3) = - 
P(n(5-_i,^-E(5-_i,^))<r) 



/3t /" dfi{x) 



2 / u — x 



T 



(215) 

(216) 
(217) 
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and for w given in ()48p . we will show 



P(V^(5-_i,^-E(5-_i,^))<r) 



1 



yj2-K/l3ta. 



Wq — X 

e 



(218) 
(219) 



Assuming (I216P and (j217p . we find that 5'^_i a converges in probability to -f / and (fT98]) 

is proved. Assuming (j218p and (j219p . we find that S^.^ + / converges in probability 
to 0, and ([20^ is proved. 

To prove ()216p . we denote for x <c the function 



gp{x; w) = n ar^ 



{w - xf/^ 



(220) 



gi3{x; w) depends on n, but we suppress that dependence to economize on notation. Let w be given 
by ()46p . uniformly for all x < c 



gp{x;w) = h 0{n ^). 

2 u — X 



(221) 



Define the (n — l)-variable function 

Pn-i,i3ixi, ■ ■ ■ ,x„_i;n,w)e'"^J=i ^^(^^j;-^) maxi<j<„_i < c, 

(222) 



Gr,/3(xi,...,x„_i;u,'«;;c) := 







otherwise. 



We have (comparing with ()194p ) 



rnS™ , „\ ^n-l,l3 {Gr.,l3ixi, . . . ,Xn-i;U,W;c)) 



nE(5-_i,^) = |-E(e™^"™-.^) 



(223) 
(224) 



r=0 



For any r G M, analogous to ()206p we have 
log E„_i^/3 (Gr_^(xi, . . . , ti, u;; c)) 



= (n-l) y r5/3(x;t(;) +p(x;ti,'«;)(i/i(x) + - 1 ) / r5r/3(x; ui) + p(x; u')di/(x) 
2 

+ ^^(^9/3(2;; w^) + p{x; u, w)) + 0(1) 



(225) 



■■logZn-i,p{u,w;c) +r 
'2 



(n - 1) / gp{x;w)dix{x) 



y gp{x;w)dv{x) + ^(^^(x; w),p(x;u,tf)). 



2r^ 

+ — A(g^(x;u;)) + 0(1), 
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where we use (j20ip in the last Hne. Using the Cauchy-Schwartz inequahty, we find (for notational 
simplicity, we write Gr^pixi, . . . , u, w; c) as Gr^/j if there is no confusion) 

^logE„_i_^(Gr,/3(a;i, . . . , x„_i; u, ix;; c)) 



-En-l,/3(Gr,/3) 



/n-1 



\k=l 



(226) 



/n-l 



f n— 1 



\k=l 



vfc=l 



>0. 



Hence log E„_i^^(G.r,/3(xi, . . . , Xn-i;u, w; c)) is a convex function in r. For any e > 0, by (1224p 



log En-l,l3{G-e,l3) - Zn-l,l3{u, W] c) ^ d log (G^./j) 



r=0 



< 



logEn^l,p{Ge,p) - Zn-l,l3{u,W;c) 



(227) 



Taking e ^ 0, by (|225D . ([2271) and ([221]) we have 

nE(5;^_^^) = (n-1) y ^/^(x; u;)(i/^(x)+ - 1^ j gfs{x;w)du{x) + ^{gi3{x;w),p{x;u,w))A + o{l), 

(228) 

and by (1220]) 

E(S;r_i,^) = 1 5^^; w)df,{x) + 0(1) = - Y / ^ + o{l). (229) 

To prove (|217p . we consider the moment-generating function of n(5'^_-^ ^ — E(5^_-|^ ^)). By ([223p 
we have 



M„(5»_^_^-E(5-_,,,))(e) =E(exp[en(5-_i,^ - E(5-_i,^))]) 

_ En-i,/3(Gg,/3(Ai, . . . , Xn^i;U, W] c)) 

exp[^n E(5;f_^_^)]i„_i,/3(tt, w; c) 



(230) 



Then by (1225]) and ([228]) we have 

^n(5™_i,,-E(S-_,,,))(0=exp 

= exp 



2f2 

-|-^(9/3(x;«;)) + o(l) 
-u — X 



(231) 



where in the last step we use (I22ip . The convergence of moment-generating function (I23ip implies 

dnzD. 

To prove (|218p . we denote for x < c the function 



1 1 

gp{x]w) := —=gp{x\w) = ^/naig- 

Vn (w — x)P'^ 



(232) 
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gj3{x;w) depends on n, and we suppress the dependence to economize on notation. Let w be given 
by ()48p . uniformly for all a; < c 



gp{x]w) = — ^- h 0(n ^). 

1 Wq — X 



(233) 



Like (j228]) . we have 

\/nE(S';^„i^^) = (n- 1) y gfi{x;w)dfi{x) +0-1^ j g,i{x]w)dv{x) 



+ -^{9p{x;w),p{x\u,w))a +o(1), (234) 



and by (|234p we obtain 



E(5'.;^_i,/3) = Vn gi3{x;w)dfi{x) + o(l) = — -^/n 



/3t ^ f dfj,{x) 



Wo — X 



(235) 



and complete the proof of (1218p . 

To prove (j219p . we consider the moment-generating function of \/n{S^_^ ^ — E(5^_^^)). Like 
()230p and (]23ip . we have the convergence of moment-generating function 



M^(5-_^_^_E(5»_,,,))(e) =E(exp[CV^(5-_i,^ - E(5-_i,^))]) 

= exp 



2£2 

-^A{gp{x;w))+o{l] 



(236) 



: exp 



2(^2 



which implies (I219p . 



□ 
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A Contour integral formula of the joint p.d.f.s of the eigenvalues in 
/3-external source models with = 1,2,4, and /3-external source 
model for all /3 > 

The goal of this appendix is two- fold. We prove Proposition 12.21 and also propose the definition of 
the /3-external source model. 

Remark A.l. The strategy in this appendix has appeared in [25, Appendix] independently for the 
purpose of proof of \25\ Theorem 1]. Since we are concerned with /3 = 1, 2, 4 cases and furthermore 
all /3 > 0, we give full detail in this appendix. 
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By change of variables and calculation of Jacobian (cf. [24^ Chapter 3]), it follows from ([T]), 
([2]) and ([9]) that the joint p.d.f.s of the eigenvalues Ai, . . . , A„ of M in the three /3-external source 
models (/3 = 1, 2, 4) are given by 

1 " r 

p.,^(Ai, . . . , A.) = 7;^|A(Ai, ...,Xn)\'J[ e'^^^^^^ / e"«T'^(^-''«^"«-^)<iQ, (237) 

where Cn^p is a normalization constant and C„_^/C„_^ is a constant depending only on n. The 
integral in (j237p is with respect to the Haar measure of the compact group G^(n), which is the 
orthogonal group 0{n), the unitary group U{n) and the compact symplectic group Sp{n) for 
/3 = 1,2,4 respectively. The matrix A„^^ is defined in ([6]) and ([TO]) , and 

A„ = diag(Ai,A2,...,A„). (238) 

Recall that in combinatorics, a partition n = (ki, K2, . . . ) is a sequence of non- negative integers 
in decreasing order, and containing only finitely many non-zero terms. We denote as the 
number of non-zero terms of k, and write k h A; if Yl\=i ~ ^■ 

Jack polynomials ci°^(xi, . . . ,x„) are n- variable symmetric polynomials indexed by partition 
K and the parameter a. For general references of Jack polynomials, see [23j and [31j. In this paper, 
we take the "C" -normalization of Jack polynomials [16], such that 

C(")(xi,...,x„) = (xi + --- + x„)^ (239) 

K\-k,l{K)<n 

The Jack polynomials with parameters 2, 1, ^ are Zonal spherical functions. See [23^ Chapter VII]. 
ci^^ are the well known Zonal polynomial in statistics [27], c2'^ are the complex Zonal polynomials, 
and are better known as Schur polynomials, and cj^^"^^ are the quaternionic Zonal polynomials. 
The integral in ()237p can be expanded in Jack polynomials: 

Proposition A.l. Let /3 = 1,2,4 and Gp{n) he 0{n), U{n) and Sp{n) respectively. If A^^p is 
defined by ^ and pOj) and A„ is defined by (j238p . then 



/ nSRTr(A„ aQAnQ'M^r) — ^"^ ' {.0-1 O-n )^ n [.M, ■ ■ ■ , 

jQ^G^in) ^^fc^;(^)<„ C^/''^(l,...,l) 



Proof. Any n-variable symmetric polynomial / can be regarded as a polynomial function from the 
spaces of n X n matrices M„(F) to F, where F stands for the division algebras M, C and H. For 
F = M and F = C, the definition is simple: If M G M„(R) or M G M„(C) and the eigenvalues of 
M are ^1, . . . , then [23l Pages 420 and 443] 

/(M) = /(ei,...,e„). (241) 

For F = H, the definition is more complicated and the reader is referred to \12>\ Page 452]. In all 
the three cases, identity ()239p implies that 

^Tr{An,pQKQ-^f=(^ C(2/'5)(A„QA„Q-^), (242) 
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where A„ = ^A„_^ = diag(ai, . . . , a„) as defined in dH). Furthermore, by general theory of Zonal 
spherical functions (e.g. [20, Proposition 5.5]) 

/ Ci^/^)(A„QA.Q'^) = • ' • • • ' ^"^ (243) 

After expanding en-RTr{A„,pQknQ ^) i^^^ power series of KTr(A„^/3QA„_/3(5""'^), we prove (|240p by 
(I2i2]) and (l2i3D . ' ' □ 

In case that A„ = diag(a, 0, . . . , 0), (I240p is much simplified by the property of Jack polynomials: 

Proposition A. 2. \31\ Proposition 2.5] If the number of nonzero variables among ai,...,an is 

less than 1{k), then C^\ai, . . . , a„) = for any a > 0. 

Therefore, in the case A„ = diag(o, 0, . . . , 0), 
Cf^)^\a, 0, . . . , 0) and cf^^{^\l, . . . , 1) can be calculated explicitly [H Table 5] 

A;— 1 2 • 

^(2//3)/„ n n^ _ „fc r-(2//3)n i X _ TT 



Cfff )(a,0, . . . ,0) = a^ ^{,7^(1, . . . , 1) = n T-l^- (245) 



Thus we have 



By [3H Proposition 2.1] and the conversion between the "J" -normalization and "C" -normalization 



|16t Table 6] , we have the identity of formal series in a 



fc=0 j = l 



Hence we obtain by Cauchy's integral formula and (I245P 

ni=o(l + fj) (2//3).. X ^_ 1 /tT 1 



(|)^fc! • • • ' - ^ n (i_,A,)/^/2^' ^^^^^ 

where the contour is taken to be a small circle around such that all Aj^ (j = 1, . . . , n) are in the 
exterior of the contour. By ()245p and ()248p . we obtain 



(|n)^q,y^(a,0,...,0)q,7^(Ai,...,A„) 

4r^(i,...,i) 



oo 



U%jr^M-'t{-^)^¥ ''''' 
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Note that (f2i9l) is valid for all j3 > 0. 

Suppose /3 is a positive number. Let m be an integer and ^ G (0, 1], such that 

^n = m + C, (250) 

we have 

1 / 1 / (^) °° 1 f —an\' \ 



/3, 



(fan)'" 27ri Jo y (1 - zAj- 



n 



(fan)™ 7^ I (u; - A^)^/ 

where the contour is large enough so that all \j are in its interior, and is in its interior if ^ 7^ 1. 
Here 

{c)i := c(c + l)(c + 2) . . . (c + n - 1) (252) 
is the Pochhammer symbol ("rising factorial"), and 

^(-i'-2,.):=E7S^-^ (253) 

is the Kummer's (confluent hypergeometric) function. See [T! 13.1.2]. Alternatively, for ^ 1 

M(l, z) = {^- l)z^-^e'j{^ - 1, z) (254) 

where 7(5,2) is the incomplete gamma function (cf. [TJ 6.5.12]), and for ^ = 1 

M(l,l,2)=e^ (255) 

We note that in the cases /3 = 2, 4, or in the case that /3 = 1 and n is even, = 1 and m = — 
Thus by (^71) . (1^ . I^m and we have 

n ^ .. I n ^ \ 

Pn,/3(Ai,...,A„) = a,/3|A(Ai,...,A„)|^ne-§"^(^^)— ^ 1 e^^^du; (256) 

where Cn^p is a constant, and the contour encloses all \j in its interior. 
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If a > and maxi<j<„ Xj < u, the contour in (|256p can be taken as Sf^ defined in (f5T]) or 
Hg defined in (I52p . where z > u. From the joint p.d.f. of Xj, it is straightforward to find the p.d.f. 
formula ([56|) of the largest eigenvalue ^max("')- Thus Proposition 12.21 is proved. 

If /3 = 1 and n is odd, similarly we obtain that the joint p.d.f. of The eigenvalues in n-dimensional 
1-external source model is 

Pn,i{Xi, . . . , A„) = 

■iny(A,) J_ / / TT I ) „ -1/2; 1 1 



C„,|A(A„ . . . , A„)| n ^"^"^^'-^^ £ [ n _ I -anw). (257) 



In Section [2] we compute the limiting distribution of ■^max('^) based on ()256p . Since the asymp- 
totic property of M(l, |,z) is similar to that of = M(l,l,z) for large z, we can compute the 
limiting distribution of Cmnxin) based on (|259p by the same method that we use in Section[21 Hence 
we can prove that Theorems 11.11 [L2l 11.31 and 11.41 hold when /3 = 1 and n is odd. 

Inspired by the Coulomb gas interpretation of the distribution of eigenvalues in random matrix 
models (see [IS])) we generalize the /3-external source model to any /3 > as the probability 
distribution of n points on the real line, such that 



p„,;3(Ai, . . . , A„) = ^|A(Ai, . . . , A„)|^ n e-^'^^^'^) 



Cn,f} i=i 



^ k\ ^ ^(2//3)/-| -.N ^ \ ) 

fc=0 Kl-fc,i(K)<n '-'K lJ-)---)J-J 

where 1^ is the potential and ai,...,a„, are external source parameters. By (j237p and (j240p . 
(j258p gives the distribution of eigenvalues of the random matrix models with external source with 
/3 = 1, 2, 4. But for other value of /?, it has no matrix interpretation. By Proposition ! A. 21 (12451) (I247P 
and ()25ip . we find that if one external source parameter is a and all others are 0, the distribution 
of the right-most point in the general /3-external source model is 



Pn,/3(Ai, . . . , A„) = C'„,;3|A(Ai, . . . , A^l^ n e-^"^(^^-) 



X 



where ^ is defined by ()250p . It is of interest to compare this formula with the rank 1 spiked Gaussian 
and Laguerre /3 ensembles studied in fS]. In the very recent preprint [19) . Forrester obtained similar 
formulas for /3-Wishart ensembles. 

B Computation of M-^iu) 

In this Appendix, we show that when /3 = 2, the •pj^2{p) in Theorem 11.31 are the same as the ]}f\(oL) 
in [21 Formula (52)], and the pj^2(«) in Theorem 11.41 are the same as the 'p'i\(a) in [3, Formula 
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(63)]. Hence we verify that the result obtained in this paper agrees with the result in [3]. Since 
Pj^2(a) and pj^2(a) in our paper are defined by G{u;a) and M.2{u) in the same way that pj_2(a) 
in [3 Formulas (52) and (63)] are defined by G{u;a) and A^i,„(u), (see (fWD . (fT68]l . ifTHT]) . (fTMIl 
and (jl83p in this paper, and [3, Formulas (177), (178) and (182)] and the comments in [3j below 
[31 Formula (52)]), we need only to show that J^2{u) / M.i^n{u) is a nonzero constant for u > e. (e 
is the right edge of J, the support of the equilibrium measure //.) The explicit formula of A^i,n('u) 
is given in [21 Formula (315)] in the case that the support of equilibrium measure is one interval. 
In this appendix we obtain that for u > e 

M2{u) = C{^{u)-j{u)-'), (260) 

where C is a constant independent of u, and, if J = [bi, 62], 

,(„)=(!i4lV". (261) 



^U-b2 J 

Thus we prove the statements above. 

To make the notations simpler, we assume J = [—1, 1] in the proof of (|260p . The generalization 
to arbitrary J is straightforward. 

From formula ([MD, (flMjl and (fT88]l . we have that 

-1 



M2{u) =exp 



A{-V{x) + log(ti - x)) 



-V{x) + log(n — x)dfi{x) 



1 v'(b)\/1- 
1 s — x 



-dsdx 



X e 



X e 



1 rl log(u-x) rl V\f)VnEdsdx 

i^-'-^ V^l^ s-x "^^"^ 

-3^ /-I ^^=f P-v- /-I (u'^^)^.) '^^^^-J-i login- x)df,ix) ^ 

1 rl log(n-j:) pi Vl-3^ 

f^J-i v^Trz2 



(262) 



(li — s)(s — x) 



dsdx 



The right-hand side of (|262p is divided into the product of three terms. The first one is a constant, 
and we compute the other two terms below. 

First we compute the third term in (|262p . Exchanging the order of integration, we have 

log(n — x) =^P-v- / 7 77 T '^^'^^ ~ T~2 / F{u,s)ds, (263) 



1 



4vr2 
where 



x" 



-1 



F{u,s)=p.v.J log 



{u — x)- 



zdx. 



(264) 



To evaluate F{u, s), we note (with the change of variable x = sin 9) 



d_ 

du 



F{u,s) =p. V. y 



1 


r 


s — u 








1 


TT 


s — u 


iVu2-l 



(u — a;)(s — x)\/l — 
1 



-d^ — p. V. 



2L s — sm ( 

2 



(265) 



TT 



(s — u)\/v? — 1 
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On the other hand, 

F{u,s) =logup.v. y 



-1 [s - X 



zdx + p. V. / log 



u — X 



1 



u 



(s — x)y/l — 



dx 



=0 + p.v. / lo^ 



u — X 



1 



(266) 



zdx. 



u {s - x)\/l - x"^ 
and it imphes that F{u, s) — )• as u — )• oo. Thus from (j265p and (j266p 



F(n,s) 



TT 



MS — 1 



arcsm s — arcsm ■ 



u — s 



Because 



arcsm s , vr. ,2 -.n, t ^ + — 1 



ds=—\og(u -l) + 7rlog 
I u — s 2 



1 arcsin ^^^^ 



_! li — s 
By (I263|) . (12671) . (1268)1 and ([269)1 we find 



-log(M -Ij-vrlog . 



(267) 

(268) 
(269) 



,4^2 J 



u + lV/^ fu-l\ 



u-l 



+ 



1/4 



U+1 



(270) 



Next we compute the second term in (1262p . For the equihbrium measure dfj,{x) = ^{x)dx on 
its support [—1, 1], By [15l Formula 6.135], we have 



^(x) = m+{x) = -m^{x) = -{G+{x) - G_(x)), 
where G{z) is an analytic function [151, Formula 6.141] 

^/^^ V'{s) 



G{z) 



{s-z)VT 



zds 



(271) 



(272) 



is an analytic function in C\ [—1, 1] and ^/ z"^ — 1 is analytic in C\ [—1, 1] with z"^ — 1 ^ z, z ^ oo. 
Thus by (I27ip . ()272p and exchanging the order of integration. 



j \og{u - s)dfj,{x) J log(u - s)\/l - x2 p. V. y 



V'(s) 



-dsdx 



1 T// 



V'(s) 



(273) 



--G{u, s)ds, 



where 



2^2 

G(n, s):=p.v. / log(n — x)\/l — dx. 

J s -X 



(274) 
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s) = tt 



I— U^^Ju?■ -\ 

yw^ — Y — u^s log ^ 



I su — \ I 

- V 1 ~ s arctan — ■ + v 1 ~ s arctan 



Therefore 



H / 1/ (sj arctan , — arctan — , | as 

2tt J_i \ 



su — 1 



Using (|275p and integration by parts, 



— n / p. V. / -dsax 

' G[u,s)as 



= --{V{l)-V {-!)) + — C V'{s) \ arctan , ^ = + 2 arctan J j"" ^ — 4 
Similarly, 



us — 1 
arcsm s — arcsm ■ 



47r2 J_^ \ u — s 



Thus by (|276|l . (|277D . ([278]) and the identity 



us-1 s (u + l)(l-s) 

arcsm s — arcsm h arctan , + 2 arctan \ / 

u-s Y(u-l)(l + s) 

S su — 1 TT 

— 2 arctan , + 2 arctan 



we have 



s^)iu'^ -1) 2 



1 log(n-x) V'{s)Vl-s 



47r2 7_i ^ 



2 

p.v. / — ^^-^ —dsdx 



p.v. / — -dsax— I log[u — x)dfj,[x) 



^v^r^d.-i-iog^i±^^r ^ 



s 



2 
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By the property [15^ Formulas 6.143 and 6.144] 

"1 V'{x) 



1 Vl — 



diJ.{x) = 0, (281) 



xV'\ 



X 



x^ 



d^{x) = 27r, (282) 



we further simplify the second factor on the right-hand of ([262]) as 

7^J-i ,JYZ^ P-^-J-i s-x 

■ ^'rr — P-v- < ^^T" — Tdsdx— f^, \oQ:(u—x)dn(x) I , / — , 

gl. 3^J-i ^77^12^^ J-1 J-1 SI ^ =2(u- Vw^-l). (283) 

Substituting (f270]l and ([283]) into ([262]) . we obtain 

I 1 \ 1/4 / 1 \ 1/4' 

u+l\' /n — 1^' 



M2(n) 



2/. ^P.v./^, ^^^^^ 



u— 1/ \n-|-l 

and prove (I260p in the case that J = [—1, 1]. The general case can be proved by a simple rescaling. 
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